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ABSTRACT 

We present a new method for the simultaneous calculation of the cosmic ray 
ionization rate, (yl^, and the ionization fraction, Xe, in dense molecular clouds. A 
simple network of chemical reactions dominant in the creation and destruction of 
HCNH"*" and HCO"^ is used in conjunction with observed pairs of rotational tran- 
sitions of several molecular species in order to determine the electron abundance 
and the Ha"*" abundance. The cosmic ray ionization rate is then calculated by 
taking advantage of the fact that, in dark clouds, it governs the rate of creation of 
Ha"^. We apply this technique to the case of the star-forming region DR21(0H), 
where we successfully detected the ( J = 3 ^ 2) and ( J = 4 ^ 3) rotational tran- 
sitions of HCNH"'". We also determine the C and O isotopic ratios in this source 
to be ^^C/^^C = 63 ± 4 and ^^O/^^O = 318 ± 64, which are in good agreement 
with previous measurements in other clouds. The significance of our method lies 
in the ability to determine A^(Hj^) and Xe directly from observations, and esti- 
mate accordingly. Our results, (n^ = 3.1 x lO^^^s^^ and Xc = 3.2 x 10~^, are 
consistent with recent determinations in other objects. 
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1. Introduction 

The ionization fraction (or the electron abundance), Xc = n{e)/ ri(H2), where n{e) is the 
electron number density and ri(H2) the hydrogen molecule number density, plays a key role in 
the chemistry and dynamics of dense molecular clouds. Due to the low temperatures within 
dense clouds, the chemistry is dominated by reactions between ions and neutral species. 
Furthermore, if a magnetic field permeates the cloud, the ions interact directly with the 
field through the Lorentz force, while the rate of coUisions between the ions and the neutral 
species determines the degree with which the neutral gas couples to the magnetic field. 
The ionization fraction therefore determines the strength of this coupling, and hence the 
ambipolar diffusion timescale in a molecular cloud, which is a measure of the stability of the 
cloud against gravitational collapse (Shu 1992). Moreover, determination of the fractional 
ionization can be essential in precisely estimating the magnetic strength in molecular clouds 
(Li & Houde 2008). 

In order to determine the ionization fraction in a cloud, it is important to identify the 
main ionization processes in the region under study. Star formation occurs in dense cores, 
which are regions of high extinction where self-shielding prevents the UV photo-ionization 
of H2. It is also expected that X-ray ionization is only significant in the vicinity of strong 
X-ray sources such as active galactic nuclei (McCall et al. 1999) and OB stars. Therefore, 
cosmic ray ionization is beheved to dominate photo-ionization in dense cores (McKee 1989). 

Cosmic rays also heat the interstellar gas and drive interstellar chemistry in dense molec- 
ular clouds. Direct determination of the cosmic ray ionization rate, (h^, is achievable by 
studying the abundance of H3+ due to the relative simplicity of its chemistry (McCall et al. 
1999). This important molecule is created via cosmic ray ionization of the H2 molecule, and 
is highly reactive with electrons and neutral species present in such clouds. For example, it 
reacts with HCN, CO and N2 to form HCNH"*", HCO"*" and N2H"'', respectively. The resulting 
ions react with different neutral molecules, producing other ionic or neutral species. For this 
study, we assume that H3+ is mainly produced by cosmic rays. 

Various techniques have been developed in the past to estimate the ionization fraction 
and cosmic ray ionization rate. For the prior, there have been studies that based its de- 
termination on measurements of the degree of deuterium fractionation through DCO^ and 
HCO+ abundance ratios (e.g., Caselli et al. 1998). The application of this technique in cold 
clouds is limited, since the freeze-out of molecules onto the grain surfaces affects the degree 
of deuterium fractionation independently of the ionization fraction (CaseUi et al. 1998). Ob- 
servations of H3+ absorption fines (McCall et al. 1999) together with appropriate chemical 
models have been used for the determination of iii diffuse clouds and envelopes of molec- 
ular clouds (van der Tak & van Dishoeck 2000; Indriolo et al. 2007). However, this method 
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is not applicable to cold starless cores since H3+ has no electric dipole moment; it cannot be 
used to trace such regions. 

In dense cores, typical values for are estimated to be (1 — 5) x 10"^^ s~^ (Dalgarno 
2006), while the fractional ionization is found to he within a small range, 10~^-^ i$ Xe 10~^-^ 
(Bergin et al. 1999). However, some recent studies have found lower values for these param- 
eters in dark clouds. For example, Maret & Bergin (2007) obtained a fractional ionization of 
5 X 10^*^ with respect to H nuclei corresponding to a cosmic ray ionization of (1 — 6) x 10^^^ 
s~^ in the Barnard 68 prestellar core, while Flower et al. (2007) reported ric/nu = 1 x 10^*^ 
and (11,, = 2 X 10~^^ in TMC-1. Since the two parameters are interdependent, one is 
usually determined by restricting the other through fitting theoretical models to observations 
(see for example, Wootten et al. 1979; Plume et al. 1998). 

In this work we show that it is possible to estimate the cosmic ray ionization rate and 
fractional ionization simultaneously, using observational spectroscopic data together with a 
simple network of chemical reactions that involve H3+ and the electron abundance. These 
reactions are responsible for the formation and destruction of well-studied molecular species 
that co-exist in dense molecular clouds. Accordingly, we chose to study HCNH+ and HCO+ 
and applied this technique to the star- forming region DR2 1(011). 

We explain the observational procedure in §2 and present our technique in §3. Our 
numerical calculations and results are described in §4, followed by a discussion and summary 
in §5. 

2. Observations 

2.1. Source description 

DR21(0H) {a = 20*^39^01^ and S = 42°22'37.7", J2000) is located about 3' to the 
north of DR21, a well-known massive star formation site in the Cygnus X region that lies 
at a distance of approximately 3 kpc (Genzel & Downes 1977). Also known as W75S, it 
is made up of several compact sources, namely DR21(0II)Main, DR21(0II)N, DR21(0II)S 
and DR21(0H)W, which are all active star-forming regions (Curran et al. 2005). This 
work is focused on the brightest component, DR21(0H)Main. Continuum studies have 
detected two dense cores, MMl and MM2, with a total mass of about 125 Mq in the center 
of DR21(0H)Main (Woody et al. 1989). This source has been extensively studied in the 
infrared and also mapped in CO (Dickel et al. 1978; Magnum et al. 1991), and no young 
stars with strong radiation fields have been observed around it (Davis et al. 2007). Also, no 
centimeter-wavelength continuum sources have been detected within the source, suggesting 
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a lack of HII regions (Johnston et al. 1984), implying that it is in an early stage of evolution 
and thus suitable for our analysis. 

2.2. Spectroscopic data 

We obtained all observations at the Caltech SubmiUimeter Observatory (CSO), located 
on Mauna Kea, Hawaii. In Table 1, we list the detections of different rotational transitions 
from the molecular species that form the main creation and destruction routes for HCNH+ 

and HCO"*" (see §3). Among these molecules, HCO"*", CO, HCN and HNC were optically 
thick, and therefore not suitable for our analysis. Instead, wc observed the optically thin 
isotopologucs Hi''^CO+, H^^cn and HN^^q^ 12,3180, ^^C^^O, and ^^C^^O in order 

to determine the ^^C/^^C and ^^O/^^O isotopic ratios needed for the calculation of the 
abundance of the main species from their observed isotopologucs. 

In October 2006, we detected HCNH+ in the J = 3 ^ 2 and J = 4 ^ 3 transitions in 
DR21(0H)Main using the 200-300 GHz receiver. For these observations, standard telescope 
efficiencies of 66 % at 222 GHz (beam width -33") and 60 % at 296 GHz (beam width ~ 
25") were used. Data were taken in a position switching mode and the pointing was checked 
regularly using scans on Uranus. H^^CO+ and H^^CN were observed during the months of 
October, November and December 1999, using the 200-300 GHz and 300-400 GHz receivers. 
These spectra were calibrated using scans made on planets available during that period 
(Mars, Jupiter and Saturn). The telescope efficiencies were ~ 70 % for the 200-300 GHz 
receiver (beam width of ~ 32") and ~ 60 % for the 300-400 GHz receiver (beam width of 
~ 20") (Houde et al. 2000). 

We obtained the final set of observational data required for this analysis in October 
and November 2007. During that period, the ^^G^^O, ^^G^^O, and ^^G^^O molecular species 
were detected in two rotational transitions (J = 2 ^ 1 and J = 3 ^ 2) towards the center 
position of DR21(0H). We also attempted to detect H3O+ in J = 3o ^ 2o and J = 32 — ^ 22 
transitions at 396 GHz and 364 GHz, respectively. Although we were successful in detecting 
the 364 GHz transition, we did not record an acceptable detection for the 396 GHz transition. 
Standard telescope efficiencies were also used for these sets of data (i.e., ~ 66 % at 200-300 
GHz with a beam width of ~ 33", and ~ 58 % at 300-400 GHz with a beam width of ~ 22"). 
Figure 1 shows the spectra for HCNH+ and H3O+. The spectra for the J = 3 — > 2 and 
J = 4 ^ 3 transitions of H^^GO+, H^^CN and HN^^C are shown in Figure 2, while the 
detections of the GO isotopologucs are presented in Figure 3. All the data reduction was 
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carried out using the GILDAS^ and CASSIS^ software packages. 

DR21(0H) is well known to have a multiply peaked structure, suggesting that the cloud 
is in an early stage of massive star formation (Richardson et al. 1994). Therefore, it is 
important to determine whether or not the emission from the molecules under study arises 
from the same volume of gas along the line of sight. In order to confirm the coexistence 
of the molecular species along each line of sight, it is useful to obtain velocity maps and 
analyze them together with line profiles. However, velocity mapping, especially for the faint 
lines of HCNH+ and H3O+, required a longer integration time than was available to us. 
Nevertheless, consideration of the spectra in Figures 1, 2 and 3 shows an alignment of the 
peak velocities of the detected lines near -3 km s^^. Therefore, to a good approximation, 
the observed molecular species are coexistent in the region along our line of sight. Notable 
exceptions are the J = 4 — > 3 transition of HN^^C and J = 3 — >■ 2 transition of ^^C^^O in 
Figure 2 and 3, respectively. The corresponding observations are, therefore, not used in our 
analysis. 



3. Method 

In this section, we describe the main creation and destruction paths of HCNH+ and 
HCO+ in dense clouds and the method used to calculate n(e), n(H^) and Ch2- In dense 
clouds, the chemistry is dominated by ion-neutral reactions and the main formation reactions 
for HCNH+ are as follows (Schilke et al. 1991) 



H+ + HCN(HNC) — > HCNH+ + H2 (1) 
HCO+ + HCN (HNC) — > HCNH+ + CO (2) 
H3O+ + HCN (HNC) — > HCNH+ + H2O. (3) 

Likewise, in dark clouds, HCO"'" is mainly formed by the reaction of 113"*" with CO (Herbst 
& Klemperer 1973; Watson 1974) 

H+ + CO — > HCO+ + H2. (4) 



lURL: http://iram.fr/IRAMFR/GILDAS/ 
^URL: http://cassis.cesr.fr/ 
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Both HCNH+ and HCO+ are mainly removed through dissociative recombination with elec- 
trons (Schilke et al. 1991; Plume et al. 1998) 



HCNH+ + e- — ^ HCN + H (5) 

— > HNC + H 

— > CN + 2H 

HCO+ + e- — > CO + H (6) 

The rate coefficients of the above reactions are taken from the UMIST database (Woodall 
et al. 2007) and are listed in Table 2. The kinetic temperature in dark molecular clouds is 
usually between 10 to 40 K; we follow Wilson & Mauersberger (1990) in adopting a kinetic 
temperature of 20 K for DR21(0H) to evaluate these rate coefficients. 

There are several other reactions that can contribute to the creation and destruction of 
HCNH+ and HCO+ in dense clouds, but they can safely be neglected due to their low rate 
coefficients and/or the relatively low abundance of the molecules involved (e.g., reactions 
involving C2H2+, H2O+, H2CO+, HNO+, C2+, and H2S+). 

We can use the aforementioned reactions to equate the rates of formation and destruction 
of HCNH+ (equations (l)-(3), and (5)) and HCO+ (equations (4) and (6)) to obtain the 
following expressions for the electron and H3+ abundances 

_ [n(HCN) + n(HNC)][n(H+)fcj + n(HCQ+)A;g + n(H30+)fc3] 

n(HCNH+)/t5 ' ^' 

n(HC0+)n(e)A;6 

Since we have obtained all the observations necessary to determine the abundance of the 
HCN, HNC, HCO+, H3O+, HCNH+, and CO molecular species, we are left with a set of two 
equations and two unknowns (i.e., n(H3+) and n(e)). It will therefore be straightforward to 
simultaneously determine n(H3"'") and n(e) using equations (7) and (8) and our spectroscopic 
data. 

The cosmic ray ionization rate can be obtained directly from the abundance of H3+, 
which forms as follows (Solomon & Werner 1971; Bowers et al. 1969) 



H2 + CR — > H+ + e- (9) 
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H+ + H2 H+ + H. (10) 

The reaction between and H2 proceeds very rapidly (Solomon & Werner 1971) and 
is limited by the abundance of H^. The rate of formation of Ha"*" is therefore governed by 
the ionization of H2 by cosmic rays. 

The main destruction path for in dense cores is through reactions with CO (McCall 
et al. 1999), due to the high reaction rate and the fact that the latter has the highest fractional 
abundance with respect to H2 in molecular clouds. Furthermore, Ha"*" is also destroyed by 
dissociative recombination with electrons through 



H+ + e- H2 + H (11) 

— ^ H + H + H. 

Although the overall rate for this reaction is an order of magnitude lower than that 
involving CO (see equation (4)), we still include it in what follows for completeness. Again, 
assuming statistical equilibrium between the formation and destruction rate of Ha"*", we 
obtain 



n(H3+)n(CO)A;4 + n(H3+)n(e)fcn 

The main advantage of our method is that it enables us to calculate Xe and simul- 
taneously through equations (7), (8) and (12), without the need to consider more extensive 
chemical networks for each molecular species involved. As mentioned before, this is because 
the required abundances are determined directly from observations. This feature differenti- 
ates ours from previous techniques. 



4. Results 



The expressions for Xc and (u,^ in the previous section are couched in terms of abundance 
but could equivalently be written in terms of column density, which is a more natural quan- 
tity to use when analyzing observation. We calculated the excitation temperature of every 
observed molecular species (or at least for one of its isotopologues) that appears in equations 
(7) and (8) in order to calculate its column density. The C and O isotopic ratios in the cloud 
were determined with observations of the ^^C^^O, ^^C^^O and ^^C^^O molecular species, and 
were used to evaluate the column density of the common molecular species (more below) . 
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4.1. Excitation temperatures and column densities 

We determined the excitation temperatures and column densities of the observed molecules 
using a LTE approach. For this purpose, we detected two transitions for each molecular 
species and obtained the ratio of their integrated intensities (Blake et al. 1987; Emerson 
1996). The brightness temperature of a source can be written as (Kutner & Ulich 1981) 



W = To 



1 1 

f>(To/Tex) — 1 (>(To/TcMs) — 1 



-t{u) 



(13) 



where T^x is the excitation temperature, Tcmb the cosmic microwave background brightness 
temperature (2.7 K), t{u) the optical depth, and Tq — hv/k, with v the frequency. The term 
involving the background brightness temperature is insignificant compared to the other term 
and is therefore neglected. The mean optical depth can be written as a function of column 
density as 



tAt; = j T{v)dv (14) 



where Av is the extent of the spectral line, A^i is the Einstein coefficient for spontaneous 
emission (listed in Table 1) and is the column density of the upper level of a transition, 
which, in turn, can be written as a function of the integrated brightness temperature and 
mean optical depth of the line as 



h(?A. 



ul 



where the term J Tf^dv is the integrated intensity of the spectral line profile and the term 
t/(1 — e"^) (hereafter abbreviated by (3) is the optical depth correction factor (Langer & 
Penzias 1993). For optically thin fines, where r ^ 1, /3 approaches unity. Within the LTE 
approximation, the population of the levels is assumed to be thermalized and the total column 
density of the molecular species along the observer's line of sight, Ntot, can be expressed as 



hcfA^ig^ J 



(16) 



where U{Tex) is the partition function at the excitation temperature Tex, while gu and 
are the upper state degeneracy and energy, respectively. 
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In order to obtain the excitation temperature of a molecular species that exhibits opti- 
cally thin line profiles, we obtained the ratio of the integrated intensities of its two observed 
transitions using equation (16). For such molecules, we adopted l3 = 1 for the corresponding 
spectral lines and solved for Tg^, assuming it to be the same for both transitions. The optical 
depths of the lines were subsequently calculated using equation (14). For example, in the 
case of H^^CN we obtained a Tg^ of 13 K and opacities of 0.09 and 0.06 for the J = 3 ^ 2 and 
J = 4 — > 3 transitions, respectively. The optical depths for all transitions thus calculated are 
hsted in Table 1. This method could not be used for ^^C^^O and ^^C^^O, as their transitions 
are not optically thin. 

Inspection of the ^^C^^O spectra in Figure 3 reveals that the line profiles are slightly 
saturated. Furthermore, the peak velocity of the J — 3 ^ 2 transition is clearly blue-shifted 
from the J — 2 ^ 1 transition and the vast majority of the other observed lines rendering it 
unsuitable for this analysis, since it possibly arises from a different region. Therefore, care 
had to be taken toward the determination of the excitation temperature and optical depth 
of this molecule. Since the line profile for the J — 2 ^ 1 transition is relatively symmetric 
and of appropriate shape, we fitted it with a Gaussian profile. We then used the equation for 
a Gaussian line from Vastel (2006) and Newton's method to numerically solve the equation 

for Nu and subsequently for t{Nu) for a very narrow velocity interval 6v centered on v. 
The dependency of the mean optical depth on Ny, follows from equation (14). This way, 
we obtained the optical depth in different parts of the line profile to model the spectrum 
using a range of excitation temperatures (10 < Tex < 25). Every optical depth value was 
converted to a brightness temperature using equation (13) and consequently, assuming a 
beam filling factor of 1, to an antenna temperature through = r]Tb, where r] is the 
telescope efficiency at the appropriate frequency. The excitation temperature of ^^C^^O was 
taken to be the temperature of the best fit to the J = 2 ^ 1 spectrum. The line profiles 
of H3O"'" (J = 32 ^ 22) and HN^^C (J = 3 — > 2) were modelled in a similar fashion to 
obtain their excitation temperatures, since these were the only credible transitions detected 
for these molecules. 

The ^^C^^O lines are optically thick, as is made clear from their self-absorbed profiles 
(see Figure 3). In order to obtain the excitation temperature for ^^C^^O, we took the ratio 
of the integrated intensities in the wings of the hne profiles, where the line optical depth is 
small (/3 ~ 1) and applied the method described above. 

Finally, the column densities of all molecules, except ^"^C^^O and ^^C^^O, were calculated 
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using equation (16) and the values of r and Te^ as explained above. These results are 
summarized in Table 3 (the cases of ^^C^^O and ^^C^^O are considered in the next section). 
The excitation temperatures obtained for the different species fall in a narrow range of 12 
to 22 K, which is consistent with the assumption of LTE and the kinetic temperature of 20 
K discussed in 53. 



4.2. C and O isotopic ratios 

The ^^C/^^C isotopic ratio was directly determined from the ratio of the column density 
of ^^C^^O to that of ^^C^^O. The observation of these relatively low-abundance species has 
the advantage of probing denser regions over the observation of more abundant isotopologues 
such as i^Ci^O and ^^C^^O. 

For the carbon isotopic ratio, we obtain 

which is consistent with the results of Langer & Penzias (1993) in other clouds. 

In order to determine the ^^O/^^O isotopic ratio, it was necessary to compare the abun- 
dance of ^^C^^O with that of ^^C^^O. Since the ^^C^^O profiles show self-absorption, we took 
the intensity ratio, R, of the J = 2 — > 1 transitions of ^^C^^O and ^^C^^O in a common 
velocity interval in the wings, where the line opacities are relatively small (we included the 
optical depth correction, /3, for ^^C^^O in this ratio). There is an uncertainty in the obtained 
value for the column density of ^^C^^O, since R may vary over the line profile. Nevertheless, 
we examined this ratio in two different parts of each wing in order to check its consistency. 
The ^^0/^*0 isotopic ratio is determined to be 

16q /i2nx 

I^^/5^(t^)=318±64, (19) 

which is in agreement with previous studies in other clouds (Penzias 1981; Polehampton et 
al. 2005). 

The column densities of the species listed in Table 3 were then calculated from their 

isotopologues using the obtained C and O isotopic ratios. We determined the column density 
of ^^C^^O and ^^C^^O and estimated the column density of molecular hydrogen in DR21(0H) 
to be (Stabler & PaUa 2004) 
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A^(H2) ~ 3.8 X 10^N{^^CO) ~ 1.1 x 10^^ cm-^. (20) 

4.3. Fractional ionization and cosmic ray ionization rate 

We determined the H3+ and electron abundances for DR21(0H) using the column den- 
sity equivalent of equations (7) and (8), and equation (20). For the fractional ionization and 
H3+ column density we report Xe — 3.2 x 10"^ and -A^(Hj) = 5.5 x 10^^ cm~^, respectively. 
McCall et al. (1999) detected absorption lines of Hs"*" in several dense molecular clouds, and 
estimated A^(H|^) to be (1 — 5) x 10^^ cm~^. We find our calculated value for A^(Hj|') lower 
than theirs. Our result for Xc agrees well with previous findings in other sources (Bcrgin 
et al. 1999; Flower ct al. 2007). Moreover, wc estimate the cosmic ray ionization rate for 
DR21(0H)Main to be = 3.1 x 10~^^s^^ through equation (12), which agrees with recent 
results for other clouds (e.g.. Flower et al. 2007). 

It should be noted that the rate coefficients for the chemical reactions are either calcu- 
lated or measured in laboratories, and have uncertainties ranging from 25% up to a factor of 
2 (Woodall et al. 2007). The OSU database has historically been used for chemical modelling 
of dark clouds; however, the latest release of the UMIST database contains dipole reactions 
with fits specific for low temperatures and is well-suited to this type of modelling. The dif- 
ference in the rate coefficients between the two databases are within a factor of two and, in 
order to investigate the effect of this, we repeated our calculations using the rate coefficients 
from the OSU database. We find an increase in the ionization fraction and the cosmic rate 
ionization rate by a factor of two, which is commensurate with the uncertainties in the rate 
coefficients mentioned above. 

The cloud size along the line of sight is needed for the calculation of the molecular 
number densities, and can be approximately estimated through the inspection of extended 
emission maps of molecules in the cloud. But the uncertainty in the distance of the source 
adds to the inaccuracy in the determination of the cloud size. We estimated the cloud size 
from the HCN ( J = 4 — >• 3) map by Richardson et al. (1986) to be 1.7 pc. Finally, the column 
density of molecular hydrogen in dense clouds is usually calculated using standard ratios of 
CO/H2 or through extinction measurements. All the uncertainties in the aforementioned 
parameters cause to be uncertain by a factor of a few. 
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5. Discussion and Sumniciry 

We have presented a simple yet novel method to simultaneously estimate the ionization 
fraction and cosmic ray ionization rate of H2 in dense molecular clouds. These parameters 
were determined in DR21(0H) using ancillary observational data pertaining to the forma- 
tion and destruction reactions of HCNH+ and HCO+, which involve the H3+ and electron 
abundance. We have made a number of assumptions regarding this method, which will be 
discussed here. 

a) The primary assumption in this method is that cosmic rays are the dominant means 
of ionization in DR21(0H). As mentioned earlier, there has been no detection of HII regions 
in DR21(0H), which imphes there are no actively photo-ionizing stars in the cloud. Gibb et 
al. (2005) mapped the DR21/DR21(OH) region at 3 mm, 850 iim and 450 iim and identified 
several submillimeter sources within the region, which only suggest the existence of deeply 
embedded young protostars in the source. We note that the value we obtain for the ionization 
fraction is consistent with our assumption. 

b) We selected the most significant formation and destruction routes for HCNH"*" and 
HCO"'" by identifying the reactions that had the largest rate coefficients and involved sig- 
nificantly abundant molecules. There are other ionic and neutral species such as C0+, C+, 
HCN+, and H2O that take part in the HCNH"*" and HCO+ chemistry. These species are 
abundant in photon dominated regions (Savage & Ziurys 2004), but can safely be neglected 
here given the lack of photo-ionizing sources in DR21(0H), as noted above. 

c) For our calculations, we have assumed chemical equilibrium within the region under 
study. Lintott & Rawlings (2006) argue that in a rapidly evolving cloud, the dynamical 
timescale may be shorter than the chemical timescale, and hence chemical equilibrium will 
not be estabhshed until the cloud has reached a quiescent stage, i.e., passed its initial collapse 
phase. They mention this timescale to be on the order of 10^ years. If chemical equilibrium 
is not met in DR21(0H), then will be overestimated by a factor of 2 — 3 (Lintott & 
Rawlings 2006). This uncertainty is comparable with the uncertainties in the rate coeffi- 
cients, observational calibrations, and measurements of cloud size and molecular hydrogen 
abundance. 

d) DR21(0H) is known to harbor water and methanol masers (associated with the 
MMl and MM2 continuum sources), being an indication of the presence of outflows within 
the source. Magnum et al. (1991, 1992) observed this source using the VLA and found an 
excitation temperature above 80 K within a region of ~ 10" centered on MMl (VLSR ~ 
-4.1 km s~^). Gas associated with such regions would understandably have different physical 
characteristics compared to the more extended, relatively quiescent regions probed with our 
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observations. It is the latter that dominates the region we observed with our larger telescope 
beam (~ 20" to ~ 30"). Nevertheless, previous studies have established that this star-forming 
region possesses a complex structure owing to the presence of deeply embedded sources, and 
the detection of molecular masers further supports this conclusion. One could thus question 
our assignment of a single kinetic temperature of 20 K for DR21(0H). But since we obtained 
a narrow range of temperatures (12 - 22 K) for all the molecular species we observed, we 
believe that our aforementioned chosen value is justified for the spatial resolution attained 
with our observations. Incidentally, this narrow range of temperatures for all species is also 
consistent with the use of the LTE approximation for this study. 

Provided the above assumptions are met, our method will be very useful for the calcula- 
tion of Xe and ill dense molecular clouds without the need for detailed chemical models. 
However, the application of this technique can be limited due to the likely difficulty in de- 
tecting HCNH+. For example, we attempted to observe HCNH+ in other molecular clouds 
such as W3(0H), AFGL2591 and AFGL490 but could not record any credible detection. For 
these situations, other molecules such as HCS+, whose formation and destruction follow a 
similar chemistry and could be coexistent with HCO"*" may be taken into consideration. 

Moreover, it is not straightforward to precisely determine the density for which our 
values of the ionization fraction and cosmic ray ionization rate apply when one considers 
that our array of observations involves molecules with a vast range of critical densities (e.g., 
from 10'^ cm^^ for CO to ~ 10*^ cm"'^ for HCN). We had to assume that all observed species 
arc coexistent within DR21(0H) in order to carry our program; however, this is probably 
not strictly accurate. This is especially the case for CO and its isotopologues, because of 
their lower critical densities, and this brings a further uncertainty in our determination of Xe 
and ■ Unfortunately, it is difficult to evaluate the magnitude of this bias with our existing 
data. On the other hand, the good alignment of spectral peak intensities and the relative 
agreement on excitation temperatures is consistent with our coexistence approximation. 

Although our method does not take account of gas-grain interactions or deuterated 
species, which may be important for clouds at very low temperatures (e.g. Flower et al. 2007), 
its simplicity and reliance on observations render it a powerful tool for the simultaneous 
calculation of Xc and ^h2- Moreover, this method enables us to indirectly estimate the column 
density of H3+ in dense molecular clouds. This is important, since the low abundance of 
Hs"*" makes its absorption lines difficult to observe, and it is not at all detectable in the 
submillimeter regime. Furthermore, the obtained values for Xc and arc applicable in the 
calculations of ambipolar diffusion timescale and the magnitude of the mean magnetic field 
in the cloud, which will be addressed in a future paper. 
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Table 1. Data for the observed molecular transitions 



Molecular transition 


Frequency 








^ peak 


jTdV 
K km 


Velocity range 


T 




GHz 


s-i 




K 


km 


km s"^ 




HCNH+ (3-2) 


222.329 


4.61x10" 


-6 


0.062 


-3.13 


0.32 ± 0.02 


(-6.1, 0.3) 


0.003 


HCNH+ (4-3) 


296.433 


1.13x10" 


-5 


0.064 


-3.27 


0.29 ± 0.02 


(-6.1, -1.3) 


0.004 


H3O+ (32-22) 


364.797 


2.79x10" 


-4 


0.16 


-3.16 


1.38 ± 0.08 


(-6.2, 1.2) 


0.02 




261.263 


6.48x10" 


-4 


0.38 


-3.07 


2.99 ± 0.03 


(-6.7, 1.1) 


0.06 


HN^^C (4-3)^ 


348.340 


1.59x10" 


-3 


0.66 


-3.81 


6.14 ± 0.08 


(-10.0, 1.0) 




Hi^CN (3-2) 


259.011 


7.72x10" 


-4 


0.96 


-2.84 


8.67 ± 0.03 


(-9.4, 3.8) 


0.09 


H"CN (4-3) 


345.339 


1.9x10-- 


i 


0.39 


-2.70 


4.17 ± 0.02 


(-7.6, 3.5) 


0.06 


H13C0+ (3-2) 


260.255 


1.34x10" 


-3 


1.02 


-2.94 


8.37 ± 0.03 


(-7.4, 2.1) 


0.09 


H"CO+ (4-3) 


346.998 


3.29x10" 


-3 


0.56 


-2.70 


4.89 ± 0.02 


(-7.2, 1.8) 


0.07 


i3ci*^0 (2-1) 


209.419 


5.23x10" 


-6 


0.20 


-3.00 


0.91 ± 0.06 


(-6.1, 0.0) 


0.01 


"CiSQ (3-2) 


314.119 


1.89x10" 


-6 


0.22 


-3.14 


0.90 ± 0.03 


(-5.9, -0.9) 


0.01 


i^Ci^O (2-1) 


219.560 


6.01x10" 


-7 


6.60 


-3.14 


41.44 ± 0.07 


(-7.0, 1.7) 


2.19^ 


i^CiSQ (3-2)^ 


329.330 


2.17x10" 


-6 


6.06 


-3.70 


40.28 ± 0.12 


(-7.0, 0.6) 




(2-l)b 


220.398 


6.07x10" 


-7 






92.70 ± 0.14 


(-8.2, 2.8) 




"Ci^O (3-2)b 


330.587 


2.19x10" 


-6 






40.45 ± 0.32 


(-7.8, 1.9) 





"^Not used in the analysis due to shifted Vpeak • 

''Two velocity components, shows strong self-absorption. 

'^optical depth in the center of the line. 
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Table 2. Rate coefficients'' of the reactions for H3+, HCNH+ and HCO 



Index'' (i) 


reaction 


kj (cm^ s ^) 




1 


H3+ + HCN(HNC) — > HCNH+ + H2 


8.1xlO-9(T/300 K)- 


-0.5 


2 


HCO+ + HCN(HNC) — > HCNH+ + CO 


3.1xlO-^(T/300 K)- 


-0.5 


3 


H3O+ + HCN(HNC) — > HCNH+ + H2O 


4.0xlO-''(T/300 K)- 


-0.5 


4 


H3+ + CO — > HCO+ + H2 


1.7x10-9 




5== 


HCNH+ + e" — > (HCN + H) or (HNC + H) or (CN - 


hH + H) 2.8xlO-'^(T/300 K)- 


-0.65 


6 


HCO+ + e- ^ CO + H 


2.4xlO-^(T/300 K)- 


-0.69 


10 


H2+ + H2 — ^ H3+ + H 


2.1x10-9 






H3+ + e- — . (H2 + H) or (H + H + H) 


6.7xlO-«(T/300 K)- 


-0.52 



*A11 the rate coefficients are taken from the UMIST database. 

''The index follows the equation number where the corresponding reaction first appears in the text. 
'^The overall rate coefficient of this reaction is the sum of the rates of all the given branch reactions. 
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Table 3. Excitation temperatures and column densities 



Molecular Species 


Tex(K) 


Ntot(cm-2) 


HCNH+ 


22 


3.2x1013 


H3O+ 


16 


1.3x1016 




12 


4.4x10^2 




13 


9.8x1012 


H13C0+ 


15 


4.8x1012 


13^180 


20 


g.oxioi-i 


12^180 


16 


5.7x10^6 


13(^160 


18 


2.9x101^ 


12^160 




1.8x1019 
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2 



m 



- 



-I 1 1 1 1 1 T" 



-i 1 r- 



HCNH"^(4-3) xlO 



i 




J 





% HCNH''(3-2) xlO 



J 





H30^(32-22) x4 



I H30-'(3o-2o) x2 







i 



-40 



-20 20 

Velocity (km/s) 



40 



Fig. 1. — The spectra for the HCNH+ transitions J = A 3 and J = 3 ^ 2, and para H3O"'" 
(J = 82 — > 22) and ortho HaO"*" (J = 3o — > 2o) towards the center position of DR21(0H). 
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20 I 1 1 1 1 1 1 1 1 — ! — I 1 1 1 1 1 1 1 r 




-20 -10 10 20 

Velocity (km/s) 

Fig. 3.— The spectra for the J = 3 ^ 2 and J = 2 ^ 1 transitions of ^^C^^O, ^^C^^O, and 
^^C^^O towards the center position of DR21(0H). 



